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Abstract: In this note we resolve the problem of getting a state-space realization compatible 
with the deterministic state-space realization and the filtering problem. 



1. 



Problem 1. Given e M"-^"- 
and C G M"-x"-, such that 



Co e Rix"-, A e 



X 
J3 



(Aq,Co) and {A,C) are observable, 
Aq = T^^AT, 
Co = CT, 
find the similarity matrix T. 

Rewrite (2) in form 

TAo = AT 

and after vectorization, we obtain: 

Vec (TAo) = (A^ ® ) Vec (T) , 
Vec {AT) = (/„^ «) A) Vec (T) . 
Equalizing (5) to (6): 



(1) 
(2) 
(3) 



(4) 

(5) 
(6) 



(A J ) Vec (T) = (/„^ «) A) Vec (T) ^ (7) 
^ {A^ (E) - ® A) Vec (T) = 
=^ Vec (T) e ker {A^ (g> - (g> A) . (8) 

Lemma 1. Mi, M2 G M"^" are two similar matrices, i.e., 
M2 = T-^MiT, T e M"^". (9) 

Then 



rank (/„^ « M2 - Aff x /„ J = n^,^ - 4^ (10) 
<^ dim {ker (/„^ ® M2 - A/f x /„ J } n^. (II) 

Proof. Consider Td G M"==x"== a singular matrix such that 

"Ai ••• 



T^^MlTd = A 



Ao 



(12) 



• • • A J 

where A^ are the blocks of a Jordan canonical form of 

Ml. Define M (/„^ « M2 - Mf (g) /„ J £ 

Multiply M by T^^®!^^ on the left side and by Td®/„^ on 



the right side. As rank (T^ ^ ® In^ ) = rank {Td ® In^ ) = 
n-J- then rank (Af) = rank { {T^'^ ® /„J M {Td ® /„J} 
On the other side. 



(T,-i /„ J Af (Td ® /„J 
= (T^-i ® /„ J ® A//2 - A^i^ ® /„ J (Td ® ) 

- (T^-i ® ) {M^ ® ) (Td /„ J . 



(13) 



Define T^"^ := 



Til T12 ■ ■ ■ Tin 
T2I T22 ■ ■ ■ T2n 



and using the 



Kronecker product definition, hence 

Tllln^ Tl2ln^ 
T2lln^ T22ln^ 



Td^®In. = 



(14) 



and 



^ M2 



'M2 
M2 







Consequently 



(T-i®/„J(/„ 

Tllln^ Tl2ln^ 
T2lln^ T22ln^ 



•M2) = 






M2 



Afa 
M2 



(15) 



• • • Afa. 



■T11AT2 T12M2 
T21M2 T22M2 

TnlM2 T„2A'f2 



TlnM2 
T2nM2 

TnnM2. 



Ma 



and 

(T^i ® 7„ J M2) (Td ® 7„, ) 
= (T-i®M2) (Td®/„, 

= (T^^Td) M2 = M2. (16) 
Using the Kronecker product properties and T4 definition: 

= (T^^M^T^^M^) «) = A ® . (17) 

From (13), (16) and (17) and the definitions of A and the 
Kronecker product, we obtain: 

= M2 - A ® 

= diag (diag {M2)]U -^^<» ^"Oi=i • ^^^^ 
As diag {M2YjLi — h-i ® is equal to 



From the observability of (Ao, Co) , T is a unique nonsin- 

CUi 
CU2 



is of full rank and 





"Cf/r 


\ 


CU2 


v 




I 



M2 - A,/„^ 
M2 - Ai7„^ 7„^ 



M2 - Ai7„^ 

whose rank is (jix — l)ni. Hence 

rank fdiag ^diag (M2)"ii - Aj (g) 7^^^^ j = ^(na;-l)ni 

' i=l 

(19) 



fe 



and^(nx-l)ni = (nx-l)^^^ = (na,-l)nj; = nl-n^, 

because the sum of the dimension of the Jordan blocks of a 
matrix is the dimension of the matrix, i.e., Yl^^i rn = n^- 
As a result rank (M) = n| — n^- □ 



gular matrix. As a result, 

[ai a.2 ■■■ oinA = C'o 



Once T is known and since is unique, having the state 
space realization {A, B, C, D) such that the conditions of 
Problem 1 are met, it is also possible to calculate 

Bo=T-^B, (24) 

Do = D. (25) 

In this manner, we obtain another minimal state-space 
realization, {Aq, Bq, Cq, Dq) , for the same LTI. 



Back to our problem, we know that: 

Vec (T) G ker {A^ (g) 7„^ - 7„^ A) 

and 

dim {ker (A^ 7„^ - 7„^ (^A))}=nx. 
Consider Ui € M.^'^'^'^^ matrices such that 

{vec (Ui) , vec {U2) vec {UnJ} 
is a basis of ker (^q ® 7„^ — 7„^ <^ A)}. Then 
T = aiUi + a2U2 + • • • + an^Un^, 



(20) 
(21) 

(22) 

where ai € R, i = l,...,nx- To calculate i = 
1, . . . , Tlx, we only have to solve equation (22). Substituing 
T according to (22) into (3): 

Co = C (aif/i + a2?72 + ■ ■ ■ + a„^UnJ 
CUi 



= [ai a.2 ■•■ oinA 



CU2 

PUn^ 



(23) 



